Ab-initio computation of the ^"^F proton-halo state and resonances in ^4 = 17 nuclei 
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We perform coupled-cluster calculations of the energies and lifetimes of single-particle states 
around the doubly magic nucleus based on chiral nucleon-nucleon interactions at next-to-next- 
to-next-to-leading order. To incorporate effects from the scattering continuum, we solve the coupled- 
cluster equations with a Gamow-Hartree-Fock basis. Our calculations for the J'^ = 1/2^ proton-halo 
state in ^^F and the 1/2"'' state in ^^O agree well with experiment, while the calculated spin-orbit 
splitting between d^/2 a-nd ^3/2 states is too small due to the lack of three-nucleon forces. We find 
that continuum effects yield a significant amount of additional binding energy for the 1/2^ and 
3/2+ states in ^^O and ^^F. 



Introduction. Halo nuclei i.e. very fragile nuclear 
systems with a halo consisting of one or more weakly 
bound nucleons, are fascinating objects. Atomic nuclei 
with halo ground states exist at the fringes of nuclear ex- 
istence close to the drip lines. Well-known examples are 
the two-neutron halo nuclei ^He, and ^^Li, the proton- 
halo nucleus ®B, and the two-proton halo nucleus ^^Ne; 
see Ref. [2!] for a recent review. Halo states can also exist 
as excited states of nuclei with well-bound ground states. 
Halo nuclei are difficult to study experimentally due to 
their feeble nature and the often small production cross 
sections. They also provide theory with a formidable 
challenge since the proximity of the continuum introduces 
a very large number of degrees of freedom. In recent 
years, several theoretical approaches have been imple- 
mented and developed that include continuum effects and 
enable theorists to describe weakly bound states, nuclear 
halos, and unbound resonances ^-8]. 

The A = 17 neighbors around ^^O are particularly in- 
teresting and significant nuclei. First, the 3/2^ and 1/2^ 
states in ^^F are bound by only 600 keV and 105 keV, 
respectively, making the latter a proton-halo state. This 
state and astrophysically relevant reactions such as the 
"'^^F(p,7)^®Ne reaction [9] have been understood within 
the shell model embedded in the continuum Q, but an 
ab-initio description is not yet available. Second, the 
ground and excited states in ^^F and ^''O determine the 
single-particle energies of proton and neutron states with 
respect to the doubly magic nucleus ^^O, respectively. 
These energies are basic ingredients of the nuclear shell 
model, and they are also key for the understanding of 
the evolution of shell structure in the fluorine and oxy- 
gen isotopes [13] ■ Recent theoretical efforts aim at ab- 
initio shell-model calculations with a core for sd-shell 
nuclei The ab-initio computation of single-particle 



energies in ^^O and ^^F is one necessary ingredient for 
such an approach. Finally, the ab-initio approach to the 
proton-halo state in -'^''F and the 3/2+ resonances in ^^O 
and ^''F provides us with an ambitious testing ground for 
the employed method, the high-precision potentials, and 



the role of three-nucleon forces. 

In this Letter, we present an ab-initio calculation of 
low-lying states of the mirror nuclei ^^O and ^''F. The 
coupled-cluster method [12] is ideally suited for this en- 
deavor, as it is a most efficient approximation for the 
computation of ground states of doubly magic nuclei, 
and the single-particle states in odd-mass neighbors can 
be computed within the equation-of-motion techniques 
[isj. For the inclusion of continuum effects, we em- 
ploy the Berggren IJ] single-particle basis of the Gamow 
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shell model [Jl, i.e. the model space consists of bound-, 
resonant-, and continuum scattering states. Our calcu- 
lations employ the chiral nucleon-nucleon interaction at 
next-to-ncxt-to-next-to leading order (N'^LO) by Mach- 
leidt and Entcm [iB l. 

Interaction and model space. We employ the intrinsic 
nuclear Hamiltonian 



H = T-T„ 
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Here, T and Tcm denote the kinetic energy and the kinetic 
energy of the center-of-mass coordinate, respectively, and 
V denotes the chiral nucleon-nucleon interaction by En- 
tem and Machleidt ^ at N^LO. 

As some of the states we seek to compute are reso- 
nances or loosely bound halo states, we need to take into 
account continuum effects. For this purpose we use a 
Berggren representation for the proton and neutron 
Si/2, ^3/2, and dc^/2 partial waves. The Berggren rep- 
resentation is a generalization of the usual completeness 
relation to the complex energy plane, so that bound-, 
resonant-, and non-resonant continuum states are treated 
on an equal footing. The Berggren ensemble has been 
successfully used within the Gamow shell model [3 (see 
Ref. [3l for a recent review), and in ab-initio coupled- 
cluster calculations of energies and lifetimes of the helium 
isotopes ITj . In constructing the single-particle Berggren 
basis, we follow the procedure outlined in Ref. We 
diagonalize a one-body Hamiltonian with a spherical 
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Woods-Saxon potential in a spherical- wave basis defined 
on a discretized contour in the complex momentum 
plane. We employ a total of 30 Gauss-Legendre mesh 
points along the contour for each of the Si/2, c?3/2; and 
partial waves. Our converged calculations are inde- 
pendent of the choice of contour, and we checked that 30 
mesh points is sufficient to reach satisfactory converged 
results for the calculated energies and lifetimes of the 
states we consider in this work. For all other partial 
waves, the basis functions are those of the spherical har- 
monic oscillator. 

Method. The computation of the ground and excited 
states in ^^O and ^^F is a three-step procedure within 
coupled-cluster theory. First, we employ the intrinsic 
Hamiltonian ([T]) and compute the ground-state energy 
Eq of ^^O. This yields a precise reference value rela- 
tive to which the single-particle energies will be deter- 
mined. In the second step, we compute the ground- 
state energy Eq and corresponding cluster amplitudes 
for a "mass-shifted" nucleus ^^O, where the mass shift 
m m' ^ m{A + 1)/A in the intrinsic Hamiltonian ([T]) 
ensures that the correct kinetic energy of the center-of- 
mass is utilized in the third step. In the third step, we 
act with an effective one-particle creation operator (con- 
sisting of superpositions of one-particle and two-particle- 
one-hole operators) onto the mass-shifted ground state of 
^^O. This yields the energies E^ = E^ +uj^ of the states 
with spin and parity /i = 1/2"'', 3/2+, 5/2+ in the A = 17 
nucleus of interest. The difference between these energies 
and the ground-state energy of ^^O are the single-particle 
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E* -Eq. We briefly de- 



scribe the three steps in more detail. 

In coupled-cluster theory [H , 1^ 2^1 , one computes the 
similarity-transformed Hamiltonian H = e~'^ He^ for the 
closed-shell nucleus ^^O. Here, T is a sum of particle-hole 
cluster operators T = X]fc=i^fc- The fc-particle fc-hole 
(fcp-fch) cluster operator 
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is defined with respect to the Hartree-Fock reference state 
0o)- Here, and in the following, we sum over repeated 
indices. The labels i, j, k, . . . {a,b,c . . .) denote occu- 
pied (unoccupied) single-particle orbitals. The opera- 
tors Up (ap annihilate (create) a fermion in orbital p. 
In practice, we truncate the cluster expansion by set- 
ting Ta = for a > 3, and treat the triples cluster T3 
in the ACCSD(T) approximation 

cluster amplitudes tf and tfj are determined from the 
condition that the similarity-transformed Hamiltonian H 
has no Ip-lh excitations and no 2p-2h excitations, re- 
spectively, from its Hartree-Fock reference state. The 
ground-state energy is the expectation value of H in the 
Hartree-Fock reference, with small corrections due to the 
approximate inclusion of triples added. This approach is 
used for the computation of the ground-state energies Eq 



and Eq of ^^O and the "mass-shifted" ^^O, respectively. 
We employ the coupled-cluster method in an angular- 



momentum coupled scheme [22, l23|. This allows us to 
obtain well-converged results for "bare" interactions from 
chiral effective field theory (EFT) in large model spaces 
consisting of 15-20 oscillator shells. 

We wish to study the low-lying states in ^''O and ^^F. 
These nuclei differ by an additional neutron or proton 
from the doubly magic ^^O. Excited states (with domi- 
nant single-particle character) can be obtained from the 
ground state of the "mass-shifted" ^^O by action of the 
excitation operator 
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Here, it is understood that the annihilation and creation 
operators on the right-hand side of Eq. ([3]) are coupled to 
the spin and parity fi of the excited proton and neutron 
states that we seek to compute, respectively. This is 
the particle-attached equation-of-motion coupled-cluster 
method with singles- and doubles excitations (PA-EOM- 
CCSD), see e.g. Refs. [10, [111. The unknowns r"" and 
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and the excitation energies uj^ relative to the ground- 
state energy of the mass-shifted ^^O are obtained from 
solving the eigenvalue problem 



(4) 



Results. We perform a Hartree-Fock (HE) calculation 
for ^^O and obtain the reference state |0o)- In order to 
assess the role of coupling to the scattering continuum, 
we present results for the A = 17 system, starting from 
a Hartree-Fock basis derived from a Harmonic Oscillator 
basis (OHF) and a Woods-Saxon Berggren basis which 
is the Gamow-Hartree-Fock basis (GHF) [l^. For well- 
bound nuclei such as ^^O, the coupling to continuum de- 
grees of freedom is negligible. The ground-state energy 
of 1^0 differs by less than 1 keV in the OHF and GHF 
basis within both the CCSD and the ACCSD(T) approx- 
imation. We found well-converged results for the ground 
state of ^^0 in 15 major oscillator shells, and the energy 
varies by less than 0.5MeV for 26 MeV< hio < 36 MeV. 
(See Refs. | 



23| for convergence details.) At the en- 



ergy minimum fiw = 34 MeV, the ground state energy of 
1^0 is -107.6McV in the CCSD, and -120.9MeV in the 
ACCSD(T) approximation. 

Figure [T] shows our PA-EOMCCSD resuhs for the 
21|. The unknown 1/2+, 3/2+, and 5/2+ single-particle energies Egp in 'F 
as a function of fku. The data points connected by dashed 
and solid lines show the coupled-cluster results obtained 
in the OHF basis and the GHF basis, respectively. The 
horizontal lines show the experimental single-particle en- 
ergies. The underlying model space includes 17 major 
oscillator shells, in addition to 30 Woods-Saxon Berggren 
states for each of the Si/2,c?3/2i and ^5/2 partial waves. 
The results obtained in the GHF basis exhibit a very 
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weak dependence on the oscillator frequency while this 
dependence is stronger for the OHF basis. In particular, 
the energies of the ^3/2 and Si/2 in the OHF basis in- 
crease with increasing frequency of the model space. The 
states in ^^O and ^^F are well-known resonances, 
and an oscillator basis is clearly not appropriate to de- 
scribe these states. For the 5/2"*' states, we find a much 
weaker effect from the continuum. This is expected since 
the I — 2 centrifugal barrier localizes this state inside 
the barrier and reduces the coupling with the external 
scattering continuum. 
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FIG. 1: (Color online) Low-lying single-particle states in ^^F 
relative to the ^"O ground-state energy as a function of oscil- 
lator frequency hw. The data points connected by dashed and 
solid lines employ an oscillator basis (OHF) and a Berggren 
basis (GHF), respectively. The horizontal lines are experi- 
mental data. 



The coupling with the scattering continuum has a 
significant effect on the 1/2^ and 3/2"'' states of ^^F 
and ^''O. Our calculations using a GHF basis yields 
^ 1.0 MeV in additional binding energy for these states 
compared to our calculations with the OHF basis. The 
effect is particularly strong in the case of the 1/2^ proton 
halo state in ^^F, which is not even bound in the OHF ba- 
sis. Similar continuum coupling effects were found for the 
1/2"'' halo state in ^^Be Q and in the low-lying states of 
the fluorine and oxygen isotopes [1, [2^ . The lack of a cen- 
trifugal barrier and the very weak binding yield a proton- 
halo (with a root-mean-square radius of r^^s — 5.333fm 
[261] ) that is difhcult to capture in the oscillator basis. 
Our calculated binding energy for this state agrees re- 
markably well with the experimental value of 105 keV. 
This finding deserves further analysis, and we need to 
estimate the effects of the omitted three- nucleon forces. 

Within chiral EFT, the leading three-nucleon forces 
consist of a long-range two-pion exchange, a midrange 
one-pion exchange, and the short-range three-nucleon 
contact interaction 271. Three-nucleon forces are ex- 



per nucleon 23|. The effect of three-nucleon forces on 
energy differences is more subtle. Within a calculation 
based on two-nucleon forces we can, however, probe the 
effect of the three-body contact by a variation of the ul- 
traviolet cutoff A. Decreasing the cutoff employed in the 
construction of the chiral interactions renormalizes the 
two-nucleon interaction and generates short-ranged three 
nucleon forces [2^. We employ the similarity renormal- 
ization group (SRG) js^] for the generation of interac- 
tions with a cutoff A, and study the evolution of the 
excited states in ^^F as the cutoff is varied. Figure [5] 
shows that the spin-orbit splitting between the c?3/2 and 
'^5/2 orbitals increases with decreasing cutoff. However, 
the Si/2 state remains virtually unchanged as the cut- 
off is lowered to A w 3.2 fm^^. This is not unexpected 
since the structure of the dilute 1/2^ halo state is dom- 
inated by long-ranged forces, and the SRG interactions 
only change the short-range contributions. Thus, our 
result for the proton-halo state in ^^F is insensitive to 
short-range three-nucleon forces. 

Let us also comment on the center-of-mass motion. 
Ref. [2^ demonstrates that the intrinsic Hamiltonian ^ 
yields a coupled-cluster wave function that factorizes to 
a very good approximation into an intrinsic part and a 
Gaussian for the center-of-mass motion. At the cutoff 
A = 2.8 fm"-'^ we checked that this is true for the low- lying 
states in the A — 17 nuclei in a wide range of oscillator 
frequencies. 
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FIG. 2: (Color online) Single-particle energies of the Si/2, 
d-i/2, and ^5/2 states in ^^F (squares, circles, and diamonds, 
respectively) and the results for a "bare" N"^LO interaction 
(dotted, dashed, and dashed-dotted lines, respectively) as a 
function of the high-momentum cutoff A. 



Table |I] summarizes our PA-EOMCCSD results for 
the 1/2+, 3/2+, and 5/2+ states in ^^O and ^^F, and 
compares with experiment. The oscillator frequency is 
huj = 34 MeV, which corresponds to the energy mini- 
mum of the CCSD and ACCSD(T) ground-state energies 
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of "'^ O. We also show the spin-orbit sphtting between the 
1^5/2 and d'iji single-particle states. 





"0 






1/2+ 


5/2+ 


Eso 


1/2+ 


5/2+ 


Eso 


GHF 


-2.8 


-3.2 


4.3 


-0.082 


0.11 


3.7 


Exp. 


-3.272 


-4.143 


5.084 


-0.105 


-0.600 


5.000 



TABLE I: Single-particle energies of the 1/2+ and 5/2+ 
states, and the spin-orbit splitting £50(^3/2-^5/2) (in units 
of MeV) in ^''O and ^'^F calculated in a Berggren (Gamow) 
basis (GHF), and the comparison to experiment 31]. 

Let us also check the consistency of the PA-EOM- 
CCSD approximation. By normalizing the excitation 
amplitudes in Eq. ([3]) to one, we can compare the norms 
of the r°- and r"'' amplitudes and get a measure of the one- 
particle structure of the A — 17 states. For the low-lying 
3/2+, 1/2+, and 5/2+ states in ^^F, we find \r''r''\ = 0.87, 
Ir^r"] = 0.92, and Ir^r"] = 0.87, respectively. We find 
similar norms for the states in ^^O. This clearly shows 
that these states are dominated by one-particle excita- 
tions from the ^^O ground state, and the PA-EOM-CCSD 
approximation is known to perform very well in this case 

Within the GHF basis, we obtain a width for the res- 
onance states. Table HIl shows the calculated energy and 
width of the (^3/2 single-particle resonance in ^^O and 
^^F in a model space with hoj — 34 MeV. The energy of 
the da/ 2 single-particle state in ^''O compares very well 
with experiment while in -'^^F it is within 0.5 MeV. The 
calculated widths are very reasonable compared to the 
experimental values, and represent the first ab-initio cal- 
culation of resonance in an ^ = 17 nucleus. 





3/2+ 


i"F 3/2+ 




Esp 
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Esp 


F 


This work 
Experiment 


1.1 
0.942 


0.014 
0.096 


3.9 
4.399 


1.0 
1.530 



TABLE II: Computed 3/2 single-particle resonance energies 
in ^'^O and ^^F compared to data 31]. The real part Eep — 
Re[E], and the width F = 2Im[i5] are given in units of MeV. 

Conclusions. We performed ab-initio coupled-cluster 
calculations of the energy and lifetimes of the low-lying 
1/2+, 3/2+, and 5/2+ states in ^'^O and "F employ- 
ing chiral nucleon-nucleon interactions and a Berggren 
single-particle basis. The single-particle energy of the 
1/2+ proton halo state in ^^F agrees well with the exper- 
iment, and we checked by cutoff variation that this result 
is not affected by short-ranged three-nucleon forces. We 
find a reduced ^3/2-^5/2 spin-orbit splitting compared to 
experiment, and confirmed via cutoff variation that this 
is sensitive to short-ranged three-nucleon forces. The life- 
times of the 3/2+ resonances in ^^F and ^^O agree rea- 



sonably well with experimental data. Our calculations 
also show that the inclusion of continuum effects is neces- 
sary for a proper description of the studied single-particle 
states. 
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